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We study the oq and Kq light scalar mesons in 2+1 flavor full QCD. Particular attention is paid 
to fitting excited states, with an eye toward determining whether scattering states are revealed. 
An ultimate goal will be to see how dynamical quarks affect the picture outlined with an earlier 
quenched study using overlap fermions, namely, that it is the ao(1450), not the ao(980), which is 
the lowest qq isovector scalar state. 
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1. Motivation 

The light scalar mesons are not unambiguously identified in terms of quark content and their 
SU (3) / classification. There are more experimental candidates for the qq nonet than can be accom- 
modated in the simple quark model. Fig. [I] shows the current experimentally known-scalar, and 
other low-lying, mesons. 
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Figure 1: Spectrum of low-lying mesons. 



An attractive classification of light scalar mesons emerges [|J] which resolves many puzzles 
in the ordering of states and size of branching ratios. This is indicated in Fig. § : a) a tetraquark 
nonet, including the a (600), ao(980) and /o(980), below lGeV, and b) a qq nonet, including the 
<?o(1450) and 7^(1430), and an almost pure glueball, / (1710), above 1 GeV. 
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Figure 2: Proposed classification of light scalar mesons. 



Our quenched calculation with overlap fermions [Q] offered support of this view in providing 
evidence that the a (600) is a four-quark state. We also found evidence [Q] that the lowest lying 
isovector scalar qq state is the ao(1450), not the ao(980). 
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Fig. |3| shows the masses of the ao and a\ as a function of quark mass. The ao has little 
dependence on mass, especially below the strange quark mass region, and extrapolates (from quite 
close to the chiral limit with pion mass as low as 180MeV) to a mass of 1.42(13) GeV. Ghost 
would-be %v\ ' states, a quenched artifact, were identified and handled by the fit, but no evidence 
of c?o(980) was seen in this channel. Indeed it was excluded with a squared spectral weight less 
than 0.015 that of the a (1450). This may account for why the #£(1430) and a (1450) have 
approximately equal masses experimentally even though the former has a strange quark. This 
counters the conventional quark-counting rule; however, the constituent quark model may have 
limitations for light hadrons where chiral symmetry plays an important role. 
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Figure 3: ao and a\ masses are plotted as a function of m\. Also shown is the two pion mass (dashed line) 
which becomes lower than the ao around the strange quark mass region. 



The proposed picture, however intriguing, is nevertheless tentative. Early two-flavor dynami- 
cal simulations [^, |j] reported values at or above 1.5 GeV, but later two-flavor dynamical calcula- 
tions [||, [7|, ||] obtained a ground state near 1 GeV. But is this the ao or Kt]^. Dynamical simulations 
must identify both the oq and scattering state(s) before the issue is settled. 

2. Lattice Details 

We use the 2+1 flavor full QCD configurations provided by CP-PACS+JLQCD These 
have renormalization-group improved gauge action and non-perturbatively 0(a) -improved clover 
quark action. 

The lattice size is 16 x 32 with lattice spacing a = 0.12fm. The light sea quark masses have 
fc = 0.13760, 0.13800, 0.13825, for which m PS (LL) / m v (LL) = 0.71, 0.66, 0.62, with pion masses 
0.84, 0.70, and 0.61 GeV. The strange sea quark mass has K$ = 0.13760. We use valence quark 
action and masses which match the sea quark (i.e. no partial quenching). 

We compute local-local two-point correlation functions with very high statistics; for each of 
800 configurations, we use 32 different delta-function sources (requiring 32 different valence quark 
matrix inversions) well separated in space and time. 

We use the Sequential Empirical Bayes (SEB) method [JlC|], a constrained-curve fitting algo- 
rithm, to fit the ground and some excited states. 



3 



Light scalar mesons in 2+1 flavor full QCD 



Terrence Draper 



3. Effective Masses and Fits 

We look first at the heaviest quark mass for which both strange and up/down quark masses (for 
both valence and sea) are equal at K$ = K = 0.13760 (m K = 0.84GeV). 

Fig. |] (left) shows the effective mass plot for the ciq (which is also the Kq since these quark 
masses are degenerate). Fig. || (right) shows the # 2 /dof for various one-state (single-cosh) fits for 
the interval [£,16]. We get good one-state fits over plateaus as long as [5,16]. Given a starting 



interval, the automated SEB method [11] adds earlier and earlier time slices, and determines when 



new terms need be added to the fit model. 
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Figure 4: Left: Effective mass plot for the isovector scalar channel with both strange and up/down quark 
masses (for both valence and sea) equal at K = K$ = 0.13760 (m K = 0.84 GeV). All effective mass plots 
are in accord with a cosh, not exponential, model; that is, for a correlation function (with periodic boundary 
conditions) which asymptotically is saturated with a single mass term, the effective mass plot will be a 
plateau right up to the middle of the lattice. Right: % 2 /dof for one-state fits over the intervals [t ,16]. 



Final single or multi-state fits give consistent results for the ground-state mass; we obtain 
1.04(1), i.e. 1.70(2) GeV, and a first-excited state much higher. The isovector scalar interpolating 
operator can excite 7TTJ scattering states (as well as 7ZT\' and KK states). If n and r\ propagate in 
the same direction in time, then they give a contribution with energy E = +m n +E mt . At this 
quark mass, the scattering state 7tr\ has an energy which is presumably close to that of the ao (since 
2m n = 1.68 GeV). In this case, the SEB method cannot resolve the two energies with our current 
statistics. 

Next, we consider the case where the strange quark mass (for sea and valence) remains at 
K$ = 0.13760 (ni n = 0.84GeV), but the light quark mass (for sea and valence) has K = 0.13825 
{m n = 0.61 GeV). Fig. || (left) shows the effective mass plot. Fig. || (right) shows the # 2 /dof for 
various one-state (single-cosh) fits for the interval [f,16]. Again we get low % 2 /&of for a wide range 
of values of t; that is, there is a long plateau. 

Multi-state or single-state fits give consistent results for the ground-state mass of 0.83(2), i.e. 
1.35(3) GeV. The first-excited state is much higher. Using the Gell-Mann-Okubo mass formula 
3m 2 = Am 2 K — nr\. and the measured kaon mass of 0.68 GeV to estimate the mass of the r\ as 
ra^ « 0.70GeV, we expect the scattering state to have an energy near 1.31 GeV which is close 
to the measured value of the ground state energy. 
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Figure 5: Same as for Fig. 4 except that while the strange quark mass (for sea and valence) remains at 
Ks = 0.13760 (m K = 0.84 GeV), the light quark mass (for sea and valence) has K = 0.13825 (m % = 0.61 GeV). 



The dip in the effective mass plot at large Euclidean time hints of evidence of a peculiar ef- 
fect due to a behavior of scattering states on a lattice with periodic boundary conditions in time. 
If the 71 and r\ propagate in opposite directions in time, they make a contribution with energy 
Am = m ri —m n [g]. Indeed an artificially-constructed correlation function with such a "wrap- 
around" contribution can agree with the data. The dip is not statistically significant however, and 
unfortunately, further doubling the statistics did not help. We cannot confirm this effect with our 
data, but it must be kept in mind in all future simulations. Rather than regard this as another con- 
tamination to isolate and control in an effort to extract the an mass, one should take the optimistic 
view and regard this as a fascinating opportunity to measure the rj mass! 

In an attempt to better isolate the cio in the correlator by eliminating this wrap-around effect, 
we recomputed the propagators with Dirichlet (fixed) boundary conditions in time. Unfortunately, 
our statistics are poor at large Euclidean times, and we are unable to get an improved fit. 
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Figure 6: Left: Same as for Fig. 4 (left) except that while the strange quark mass (for sea and valence) 
remains at K$ = 0.13760 (m K = 0.84GeV), the light quark mass (for sea and valence) has K = 0.13800 
(m-Ti = 0.61 GeV). Right: The effective mass plot from the second trajectory of 230 configurations (green 
triangles) and that of the remaining two trajectories combined (red dots). 
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The case where the strange quark mass (for sea and valence) remains at K$ = 0.13760 (m n = 
0.84GeV) and the light quark mass is intermediate at fc = 0.13800 (m„ = 0.70GeV) is most pecu- 
liar. Fig. 6 (left) shows the effective mass plot. In contrast to the previous two cases, here there is 
no stable plateau, and the error bars are smaller than expected. We cannot make a fit in this case. 

The dynamical gauge configurations were created from three trajectories. Fig. |6| (right) shows 
the effective mass plot from the second trajectory of 230 configurations (green triangles) and that of 
the remaining two trajectories combined (red dots). The latter resembles those of the other masses, 
but the former is peculiar with relatively small error bars and lack of a plateau. This behavior of 
having one trajectory giving disparate results is not replicated in other channels. 

In summary, we are only able to obtain ground state fits in the isovector scalar channel for two 
of the three light quark masses studied. Although we don't attempt a formal chiral extrapolation 
with only two points at rather large quark mass (with m n = 0.84GeV and 0.61 GeV), the downward 
trend suggests a ground state mass near lGeV, rather than near 1.5 GeV, in this channel. At face 
value, this suggests that the lowest lying isovector scalar qq state is the £20(980), not the ao(1450). 
This agrees with that of other recent dynamical quark simulations [||, 0, ||] and at first glance 
appears to be at odds with the picture painted with a quenched overlap calculation at low quark 
mass [^]. But we cannot jump to this conclusion because the measured ground state energies are 
so close to expected Ttr\ scattering states. Dynamical simulations must identify both the ao and 
scattering state(s) before the issue is settled. 

The use of smeared sources would alter the proportion of qq and scattering states, allowing 
SEB or a variational calculation to resolve the states. We might also expect a resolution with a new 
set of 2+1 flavor dynamical gauge configurations at lighter quark masses. 



4. Hybrid Boundary Conditions 

For our quark masses, the kinematics are such that we expect there to be a 7117 scattering state 
near the measured ground state in the isovector scalar channel. To reveal this, we recomputed the 



valence quark propagators with hybrid boundary conditions (HBC) [11], using the same configu- 
rations. In our version of HBC, we use anti-periodic boundary conditions in space for the valence 
quarks, while the dynamical quarks still have periodic boundary conditions. Although the valence 
quarks are given non-zero momentum by the interpolating fields, a qq ground state will still have 
zero momentum. On the other hand, in the scattering state each meson must have non-zero mo- 
mentum (since the lowest Fock component has one valence and one sea quark). Thus, HBC raise 
the energy of the scattering states while leaving that of qq unchanged, potentially allowing SEB 
to resolve the states. (A caveat is that with mixed boundary conditions (valence versus sea), ghost 
states may complicate matters.) 

Using HBC for K$ = 0. 13760 and K = 0. 13825, we measure the energy of the lowest isovector 
scalar state to be 1.50(3)GeV. This is higher than the 1.35(3) GeV measured with periodic bound- 
ary conditions (PBC). This strongly suggests that the measured PBC ground state energy (GSE) is 
that of a scattering state, and casts doubt that the ciq mass extrapolates to a value near 1 GeV in the 
chiral limit. However, this HBC GSE does not exceed the PBC GSE by the expected momentum- 
dependent amount, if both are scattering states. Rather, the measured HBC GSE is less, and could 
be that of the qq. 
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We propose the following scenario which is consistent with our (limited) observations: at our 
lowest light quark mass the PBC GSE of 1.35(3) GeV is that of a scattering state, with the single- 
particle (ao) state lying higher. As HBC are imposed, the scattering state energy increases, leapfrog- 
ging that of the single-particle state which does not change much, leaving the single particle state 
exposed as the (new) ground state at 1.50(3) GeV. Furthermore, at higher light quark mass (degen- 
erate with the strange), we observe that the GSE does not change much from PBC (1.70(2) GeV) 
to HBC (1.69(2) GeV), suggesting that the ground state in this case is a single-particle (ao) state. 
In addition, although a naive chiral extrapolation for PBC would give an extrapolated value near 
lGeV, this extrapolation should not be done, as it mixes single-particle and scattering states! The 
HBC GSE chiral extrapolation is much less steep than for the PBC GSE, with an extrapolated value 
much higher than 1 GeV. 

Unfortunately, our evidence for such a scenario is incomplete. With the current data set, the 
statistics are not sufficient for SEB to convincingly extract excited states. We need to measure these 
to verify that we see HBC raise energies by the expected amounts. Then we can unambiguously 
identify and distinguish scattering states from single-particle states. Until this is resolved, we urge 
the community to keep an open mind about whether or not the ao(980) is a qq state. 

Scalar mesons are an increasingly rich, fascinating, and frustrating forum. 
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